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ABSTRACT 


The  author  employs  a  16K  linked,  FORTRAN  II  version  of  the 
problem  oriented  PACER  simulation  program  written  by  Prof.   Paul 
T.  Shannon,  Dartmouth  College,   Hanover,  New  Hampshire  to 
develop  a  non-linked,  FORTRAN  IV  program  which  is  compatible  with 
the  GE  625  system  installed  at  the  University  of  Kansas. 

The  author  then  formulates  a  steam  power  system  model  and 
tests  several  configurations  of  this  system  to  demonstrate 

(1)  the  efficacy  of  the  modified  PACER  program  and 

(2)  the  potential  usefulness  of  such  a  modified  PACER 
program  in  either  graduate  or  undergraduate  chemical 
engineering  curriculum. 
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Chapter  I 
INTRODUCTION 

General  .     In  early  19  65  the  National  Science  Foundation  estab- 
lished a  project  at  the  University  of  Michigan  to  evaluate  the  current 
and  future  uses  of  computers  in  engineering  design,  education.    The 
final  report  issued  by  this  project  observed  that  the  majority  of 
analytical  work  done  in  upper-level  engineering  courses  was  of  a 
routine  nature,  where  the  solution  procedures  are  generally  brief  and 
straightforward  and  where  there  is  normally  only  one  correct  answer. 
For  this  reason,  the  report  stated,   engineering  design  courses  had 
fallen  into  general  disregard  and  in  some  instances  they  have  been 
eliminated  from  the  curricula.     Rather  than  take  a  pessimistic  view  from 
these  circumstances,  the  report  concluded  that: 

Such  courses  can,  with  the  introduction  of  the  computer, 
take  on  a  new  Taster  and  make  significant  contributions 
to  the  education  of  our  engineering  graduates.     In  the 
past,  the  engineering  student  has  only  had  time  to 
investigate  perhaps  one  or  two  possible  (and  usually 
not  very  effective)  designs.    With  the  new  hardware  and 
software  systems,  the  student  designer  will  for  the  first 
time  have  the  opportunity  to  investigate  a  large  number 
of  feasible  solutions  to  his  design  problem.  .  .  . 

In  effect,  because  of  his  exposure  to  a  large  number  of 
designs  and  his  ability  to  manipulate  system  configurations 
and  design  variables  easily,  the  student  can  gain  some 
facets  of  "engineering  experience"  which  cannot  now  be 
acquired  outside  the  industrial  environment  (and  which 
would  require  a  much  longer  elapsed  time).    Thus,  in 
addition  to  its  use  as  an  analysis  tool,  the  time  shared 
computer  with  appropriate  design-oriented  software 
should  enable  the  student  to  improve  his  intuition  and 
understanding  of  the  nature  of  good  design:    this  is  the 
raison  d'etre  of  engineering  design  courses. x 

One  of  the  most  promising  approaches  to  the  demand  for  effective 
engineering  design  education  lies  in  the  recent  development  of  problem- 
oriented  approaches  to  programming  languages.    The  problem-oriented 
approach  involves  the  development  of  a  general  processor  or  executive 
program  which  is  capable  of  handling  a  general  class  of  problems  with 
a  minimal  amount  of  programming  effort  on  the  part  of  the  problem 
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designer.    This  frees  the  designer  to  concentrate  on  the  primary 
objective  of  system  design  and  analysis. 

Of  the  problem-oriented  simulation  programs  developed  thus  far, 
one  of  the  most  notable  is  PACER.     PACER  was  developed  by  Prof.   Paul 
T.  Shannon  of  Dartmouth  College  and  was  developed  primarily  to 
process  material  and  energy  balances  of  typical  system  design  problems 
encountered  in  chemical  engineering.     Solutions  are  processed  by  a 
modular  or  unit  operations  concept  with  iterative  calculations  under- 
taken when  recycle  conditions  are  encountered. 

Objective  and  Scope  of  Thesis.    This  thesis  will  examine  the 
historical  development  of  digital  simulation  programs  culminating  with 
the  development  and  some  applications  of  PACER.    The  essential 
concepts  and  characteristics  of  PACER  as  developed  by  Prof.  Shannon 
will  be  described  and  illustrated.    The  author  will  further  discuss  the 
modifications  made  to  PACER  to  make  it  compatible  with  the  computation 
system  installed  at  the  University  of  Kansas.    Additionally,  the  author 
will  trace  the  development  of  a  semi-realistic  model  of  a  steam  power 
plant  developed  to  evaluate  the  modified  PACER  program.    As  a  related 
objective,  the  usefulness  of  the  model  in  an  undergraduate  chemical 
engineering  education  curriculum  will  be  explored. 

Finally,  the  advantages  and  limitations  of  the  modified  PACER 
program  and  suggestions  for  its  further  improvement  and  refinement 
will  be  presented  for  consideration. 


Chapter  II 
HISTORICAL  DEVELOPMENT  OF  DIGITAL 
SIMULATION  PROGRAMS  AND 
APPLICATIONS  TO  ENGINEERING  EDUCATION 

Historical  Development  of  Digital  Simulation  Programs .     Initial 
attempts  toward  computer  simulation  of  chemical  engineering  systems, 
were  made  in  the  early  1950's  and  were  limited  almost  solely  to  analog 
computers.    The  digital  computer  at  that  time  required  tedious  machine- 
language  programming  and  it  was  slow  and  expensive  for  such  uses. 
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Accordingly,  it  could  not  effectively  compete  with  the  analog  computer. 

(Figure  1  presents  a  schematic  illustration  of  the  historical  development 
of  simulation  programs). 

R.  G.  Selfridge  is  generally  credited  with  the  development,  in 

1955,  of  the  first  digital  simulation  program  for  use  on  the  IBM  701 
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computer    .    This  initial  effort  was  followed  during  the  period  from  195  6 

to  1958  by  the  development  of  a  group  of  digital  computers  known  as 

digital  differential  analyzers  (D.D.A.  's) .    The  D.D.A,  's  were  designed 

to  duplicate  the  functioning  of  the  analog  computer,  but  proved  to  be 
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exceedingly  slow  and  cumbersome    . 

One  of  the  most  significant  advances  in  digital  simulation 
occurred  with  the  development  of  procedure-oriented  languages,   such 
as  FORTRAN,   in  the  late  1950's.    The  availability  of  these  languages 
facilitated  the  heretofore  cumbersome  programming  task  by  permitting 
the  programmer  to  write  his  program  in  terms  closely  resembling  the 
symbolic  representation  of  problem  variables  and  parameters.    A  com- 
piler converted  the  symbolic  program  into  machine  language  instructions. 
The  FORTRAN  language  additionally  was  provided  with  a  library  of 
functions  and  subroutines  which  performed  commonly  encountered  mathe- 
matical operations  such  as  logarithims,  exponentation,  etc.     FORTRAN 
coupled  with  increasingly  spectacular  improvements  in  memory  capa- 
bility and  speed  of  operations,  where  operations  are  currently  specified 
in  nano- seconds  (billionths  of  a  second),  gradually  began  to  assart  the 
effectiveness  of  digital  simulation  programs. 
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Figure  1.  Development  of  Simulation  trograms- 
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The  relative  advantages  and  disadvantages  of  digital  and  analog 
simulation  became  more  distinct  as  each  system  was  progressively 
developed.     It  was  readily  acknowledged  that  the  analog  computer 
enjoyed  the  advantage  of  parallel  operation  which  obviated  the  careful 
attention  to  logical  execution  sequence  and  extensive  anticipation  of 
alternatives  required  by  the  sequential  nature  of  the  digital  computer. 
Digital  computers  normally  require  extensive  reprogramming  for  different 
problem  structures.     Mathematical  integration,  a  natural  operation  for 
an  analog  computer,  requires  specific  and  extensive  programming  on  a 
digital  computer.    On  the  other  hand,  the  digital  computer  provides 
certain  specific  advantages  which  the  analog  computer  cannot  approach: 
greater  reliability,  memory  and  logic  capabilities  as  well  as  a  decidedly 
better  dynamic  range.  ° 

An  extensive  effort  was  undertaken  to  incorporate  the  best 
features  of  the  analog  and  digital  computer  into  a  digital  simulation 
capability  referred  to  as  a  digital-analog  simulator.    These  simulation 
systems  generally  consist  of  block-oriented  subroutines  which  perform 
or  simulate  the  natural  functions  of  an  analog  computer  such  as 
integration  and  differentation. 

Two  representative  examples  of  digital-analog  simulators  are  the 
MIDAS  and  MIMIC  programs.    The  most  significant  capability  of  MIDAS7 
is  its  integration  function  which  automatically  adjusts  the  integration 
step  to  a  specified  size  and  combines  this  feature  with  a  fifth-order 
Milne,  predict- correct  integration  routine.    These  characteristics 
provide  the  user  with  assurance  that  the  integration  is  being  performed 
correctly." 

The  next  effort  in  digital  simulation  involved  the  development  of 
equation-oriented,  digital-analog  simulation  programs  which  did  not 
rely  upon  the  block-oriented  structure  imposed  on  MIDAS  and  its 
contemporaries.     Chief  among  this  group  of  programs  were  DSL/90y 
and  MIMIC.10 

MIMIC,   introduced  in  19  65  by  the  authors  of  MIDAS,  differs  in 
concept  from  the  less    flexible   block-oriented   MIDAS.     MIMIC  offers 
several  other  features   not  comparably  available  in   MIDAS,   such  as 
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a  wider  variety  of  available  subroutines ,  FORTRAN  compatibility  and 
intuitive  solution  capability  for  implicit  algebraic  equations. 

A  separate  class  of  programming  languages  has  been  designed  to 
simulate  systems  characterized  by  a  sequence  of  discrete  events  such 
as  schedule  optimization  and  resources  allocation.     Representative 
examples  of  this  type  of  program  are  GASP,12    GPSS,13    SIMSCRIPT, l3 
and  CSL.13 

The  most  significant  current  advance  in  the  field  of  digital 
simulation  is  the  development  of  problem-oriented  or  application- 
oriented  programming  languages  which  are  designed  to  handle  a  specific 
class  or  specific  classes  of  problems.     Some  of  the  more  significant  of 
this  type  of  program  are: 

14  15 

1 .  BLODIB       and  ECAP  '    for  signal  analysis  and  circuit 

design  in  electrical  engineering 

2.  DYANA  "     for  vibration  analysis  in  mechanical  engineering, 

and 
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3.  PACER      '         for  system  simulation  in  chemical 

engineering. 

PACER,  an  acronym  for  Process  Assembly  Case  Evaluator  Routine, 
was  conceived  in  1961  by  Prof.   Paul  T.  Shannon  while  he  was  in  the 
employ  of  Humble  Oil  and  Refining  Company.    After  extensive  modifi- 
cation and  development,   PACER  received  its  first  practical  test  in 
October  1964  when  Prof.  Shannon  joined  with  Prof.  A.  I.  Johnson  of 
McMaster  University/   Hamilton,  Ontario,   Canada  and  Dr.  Norman 
Cooke  from  Canada  Industries  Ltd.  to  undertake  a  simulation  of  C-I-L's 
300  ton/day  contact  sulfuric  acid  plant  in  Hamilton.    This  initial  under- 
taking was  successfully  completed  and  reported  in  April  1965.    The 
sulfuric  acid  plant  simulation  involved  the  use  of  42  equipment,  infor- 
mation gathering  or  parameter  setting  subroutines  which  were  inter- 
connected by  approximately  70  information  transmittal  streams.    The 
model  thus  conceived  involves  the  solution  of  approximately  5  00 
simultaneous  equations  involving  almost  1000  stream  variables  and  2  00 
equipment  parameters.    A  complete  simulation  run  required  8  minutes  of 
computer  time  on  the  IBM  7040  computer. 
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In  1965  the  staff  and  the  senior  chemical  engineering  students 
at  McMaster  University  joined  with  C-I-L  in  the  design  of  a  600  ton/day 
sulfuric  acid  plant  based  on  metallurgical  off-gas.    The  results  of  this 
study  were  presented  to  C-I-L  and  representatives  of  other  companies 
in  April  1966.    Today,  more  than  twelve  industrial  firms  are  using 
PACER  in  the  analysis  of  their  processing  problems,  and  thirty-nine 

universities  are  integrating  PACER  into  their  teaching  and  research 

..    ...        19 
activities . 

Use  of  Computers  and  Simulation  Programs  in  Engineering 
Education.     In  1959  the  College  of  Engineering  of  the  University  of 
Michigan  undertook,  under  the  sponsorship  of  the  Ford  Foundation,  a 
study  of  the  feasibility  of  broad  scale  integration  of  electronic  computer 
use  into  the  educational  process.    The  final  report  of  this  study 
recommended  that: 

(1)  All  engineering  students  be  introduced  to  computing 
techniques  early  in  their  college  education,  preferably 
at  the  freshman  or  sophomore  level. 

(2)  That  they  develop  a  proficiency  in  a  procedure- 
oriented  language  such  as  FORTRAN,  ALGOL,   MAD,  etc., 
and 

(3)  Required  upper-level  engineering  courses  incorporate 
the  requirement  for  computer  solutions  of  several  engineering 
problems  of  graded  difficulty  by  the  student  prior  to  gradu- 
ation. 20 

A  disciplinary  appendix  to  the  final  report  entitled  "The  Use  of 
Computers  in  Chemical  Engineering  Education"  cited  the  purely 
educational  aspects  of  computers  in  the  education  process,   stressing 
the  salutary  effects  of  precise  formulation,  definition,  and  logical 
organization  required  for  computer  solution  of  engineering  problems. 

A  survey  of  engineering  school  deans  and  department  heads 
taken  in  1962  as  a  part  of  the  Project  revealed  that  approximately  30% 
of  all  engineering  students  attended  schools  which  had  a  required 
introductory  computer  course  for  undergraduate  students.    The  study 
further  indicated  45%  of  engineering  students  attended  colleges  or 
universities  which  required  students  in  some  engineering  disciplines 
to  attend  such  an  introductory  course.    A  subsequent  study  conducted 
during  the  1965-56  academic  year  showed  that  approximately  87%  of 
all  engineering  undergraduate  students  received  some  required 


3 

instruction  in  digital  computer  techniques.    An  analysis  of  this  data 

indicates  that  in  the  near  future  almost  all  undergraduate  engineering 

21 

students  will  attend  a  required  digital  computer  introductory  course. 

Quality,  however,  apparently  has  not  matched  quantity  as  the 
summary  report  on  "Computers  in  Engineering  Design  Education" 
concluded.   (See  Chapter  I).     Most  of  the  computer  solutions  assigned 
in  upper  level  engineering  courses  are  routine  in  nature.     One  major 
shortcoming  in  the  field  of  engineering  design  has  been  traditional 
reticence  of  design  instructors  to  employ  sophisticated  mathematical 
analysis  techniques  in  their  courses.     Further,  existing  computer 
hardware  and  programming  languages,  due  primarily  to  their  general 
nature,  were  not  amenable  to  flexible  use  in  design  problems. 

However,  two  recent  developments  in  the  field  of  computer 
technology  hold  the  promise  of  a  vast  improvement  in  the  area  of 
engineering  design  education.    The  first  of  these  new  developments  is 
the  time- shared  computer  with  remote  terminal  installations.    These 
facilities  will  permit  a  large  number  of  users  to  have  almost  instan- 
taneous access  to  the  computer  and  will  permit  these  users  to  communi- 
cate with  the  computer  in  an  almost  conversational  manner.    This  time- 
sharing capability,  coupled  with  the  development  of  a  new  hierarchy  of 
problem-oriented  computer  languages  and  processors  such  as  PACER, 
permits  the  visualization  of  the  investigation  of  a  large  number  of 
feasible  configurations  of  complex  engineering  problems  by  a  student 

in  the  course  of  a  semester.     In  effect,  the  computer  will  then  take 
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on  the  role  of  an  effective  imparter  of  "engineering  experience". 


Chapter  III 
CONCEPT,   ORGANIZATION  AND  OPERATIONS  OF  PACER 

Concept  of  PACER.     As  previously  mentioned  PACER  was  developed 
by  Prof.  Shannon  to  provide  a  highly- flexible,   systematic  method  for 
the  simulation  and  design  of  complex  chemical  engineering  systems  and 
processes.    The  basic  concepts  underlying  the  PACER  system  may  be 
summarized  as  follows: 

1.  PACER  is  a  systematic  method  which  permits  a  progressive 
development  of  comprehensive  models  of  chemical  engineering  systems 
or  processes . 

2.  The  processing  or  operational  sequence  and  the  physical 
configuration  of  the  system  or  process  are  included  in  the  input  data 
provided  by  the  user  and  can  be  readily  modified  to  explore  alternate 
system  or  parametric  configurations. 

3.  Mathematical  models  of  individual  equipments  are  available 
for  use  during  development  of  the  total  system  or  process. 

4.  Earlier  work  and  results  may  be  incorporated  into  the  process 
model  with  minimal  effort  on  the  part  of  the  user. 

5.  The  successful  completion  of  each  analysis  contributes  to 

the  development  of  a  library  of  mathematical  equipment  models  which 

23 

may  be  used  in  subsequent  analyses. 

PACER  conceptually  then  differs  from  a  procedure-oriented 
language  such  as  FORTRAN.     FORTRAN  permits  a  problem  solver  to 
write  solutions  for  mathematical  and  engineering  problems  in  a  language 
which  symbolically  approximates  commonly- encountered  engineering 
and  mathematical  terms.     It  does  not  provide  a  general  or  even  a  specific 
solution  method.     It's  greatest  advantage  lies  in  its  flexibility — the 
ability  to  represent  many  types  of  problems,  ranging  from  the  social 
sciences  and  law  to  nuclear  physics,  with  the  same  language. 

This  generality  also  represents  a  deficiency  of  FORTRAN.  Its 
use,  in  many  cases,  requires  the  complete  restructuring  of  complex 
solution  programs  for  only  minor  system  or  problem  configurations . 
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Such  a  requirement  can  be  exceedingly  cumbersome  to  a  systems 
designer  or  analyst  who  is  primarily  interested  in  the  results  of  the 
simulation. 

PACER  was  developed  to  provide  to  a  chemical  engineer  the 
ability  to  rapidly  and  accurately  solve  varying  configurations  of  complex 
chemical  engineering  systems.     PACER  accomplishes  this  by  providing 
a  rigidly- structured  program  which  is  capable,  with  minimal  effort  on 
the  part  of  the  user,  of  solving  complex  chemical  engineering  systems. 
No  tedious  rewriting  of  the  program  is  necessary  for  changes  in  config- 
uration.    PACER,  then,  provides  the  logical  next  step  in  computer 
language  design-- macro- procedure-oriented  programs. 

It  permits  a  user  to  systematically  describe  (by  the  use  of  a 
highly- structured  data  deck)  the  topology  and  parametric  constraints 
of  the  system  to  be  solved.     Major  changes  in  the  actual  physical 
configuration  of  a  system  can  be  accommodated  with  relatively  simple 
changes  to  the  data  deck.     It  also  provides  the  user  the  option  of 
specifying  or  not  specifying  the  order  in  which  equipment  will  be  cal- 
culated.    If  the  order  of  calculation  is  not  specified,   PACER  enters  into 
a  trial  and  error  procedure  to  determine  the  solution  order. 

The  PACER  concept  of  operations  may  be  effectively  illustrated 
by  the  presentation  of  a  simple  problem  involving  (1)  a  flow  mixer, 
(2)  a  heater  and  (3)  a  flow  splitter.     Figure  2  presents  the  flow  diagram 
for  such  a  system,  as  well  as  the  major  components  of  the  data  deck 
which  would  be  used  to  communicate  the  system  configuration  and 
parametric  representation  to  the  PACER  executive  routine.    A  more 
extensive  discussion  of  the  data  deck  and  its  components  is  presented 
subsequently  in  this  chapter  and  Appendix  A. 

Organization  of  PACER.    The  major  elements  of  PACER  are: 

1.  A  mainline  program  which  serves  as  the  general  processor 
for  PACER  operations  and  conducts  all  iterative  calculations. 

2.  Information  modifying  and  parameter  setting  subroutines. 
These  include: 

(a)     EQCAL  -  Calls  proper  equipment  subroutines  to  calculate 
system  or  process  equipments. 
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a.  Flow  Diagram 

EQPT. 

1  MIXER 

2  HEATER 

3  SPLTR 

b.  Process  Matrix 

EQPT. 
NO. 

1 
2 
3 


2 

5 

-5 


STREAM  NUMBERS 


-3 

-4 
-6 


0 

0 

0 

0 

0 

0 

0 

0 

0 

EQUIPMENT  PARAMETERS 


0 

0 

0 

0 

0 

0 

17. 

350. 

2. 

43.5 

0 

.2 

.8 

0 

0 

c.      Equipment  Parameters 

STREAM 
NO. 


STREAM  VARIABLES 


1. 

0 

15. 

80. 

100. 

0 

100. 

0 

2. 

0 

15. 

60. 

200. 

200. 

0 

0 

3. 

0 

15. 

0 

0 

0 

0 

0 

4. 

0 

15. 

0 

0 

0 

0 

0 

5. 

0 

15.    • 

0 

0 

0 

0 

0 

6. 

0 

15. 

0 

0 

0 

0 

0 

d.      Stream  Variables 


STREAM  VARIABLES  ABSOLUTE  FRACTIONAL  CONVERGENCE  CRITERIA 


1.  i.  1.  .01 

e.      Convergence  Criteria 

EQPT.       MODE 


.01 


.01 


.01 


.01 


f.      Equipment  Calculation  List 
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(b)  TEST  -  if  calculation  is  in  an  iterative  (recycle)  mode  this 
subroutine  tests  the  values  of  the  output  stream  variables  for  the  equip- 
ment being  calculated  to  determine  if  the  calculated  values  meet  con- 
vergence criteria  established  by  the  user. 

(c)  XDATA  -  provides  the  space  for  additional  equipment 
variables  if  required  by  the  user. 

(d)  MAJOR  -  called  by  PACER  when  order  of  calculation  of 
equipment  is  not  supplied  by  user.     MAJOR  directs  search  for  list  of 
equipments  which  will  calculate  as  output  streams  up  to  three  assumed 
known  input  streams.    When  such  a  list  is  developed  the  equipments  are 
calculated  as  either  direct  or  iterative  calculations  depending  on  their 
configuration. 

(e)  EQUIP  -  serves  as  a  bookkeeping  subroutine  for  all  equip- 
ments which  are  calculated  directly0    Bookkeeping  functions  for  equipment 
calculated  in  an  iterative  mode  are  performed  by  the  mainline  program. 

(f)  GUESS  1  -  forms  under  the  direction  of  MAJOR  an    ordered 
list  of  unknown  input  streams.     From  this  list,  up  to  three  unknown 
streams  are  assumed  known  and  a  search  is  made  to  find  a  list  of  equip- 
ment which  will  calculate  as  output  streams  the  assumed  known  input 
streams . 

(g)  GUESS2  -  carries  out  under  the  direction  of  GUESS1  the 
search  for  a  list  of  equipment  which  will  calculate  as  output  streams 
those  unknown  streams  assumed  known  in  GUESS1.     If  such  a  list  of 
streams  can  be  found,  GUESS2  calls  for  appropriate  calculation. 

(h)     DREAD  -  reads  data  deck  for  problem  configuration  (run). 
A  new  set  of  data  is  required  for  each  configuration  (run). 

(i)      EXREAD  -  reads  additional  equipment  parameters  if  pro- 
vided by  user  as  well  as  order  of  calculation  if  specified. 

0)      SEFLAG  -  sets  numerical  indicators  (flags)  indicating 
status  of  equipments  and  streams  provided  as  data.     Flags  may  indicate 
either  that  an  equipment  is  used  or  not  used  or  that  a  stream  is  a  feed 
stream,  product  stream,  internal  stream  or  not  used. 

(k)     DPRINT  -  prints  data  supplied  by  user.     Primarily  used  as 
a  debugging  feature. 
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(1)      EXPRNT  -  prints  additional  equipment  parameters  supplied 
or  calculated  by  PACER  as  well  as  order  of  calculation  specified  by  user. 

(m)     RPRINT  -  prints  results  of  run  conducted  by  PACER.     Results 
are  printed  at  completion  of  each  run. 

3.  Equipment  subroutines  which  are  mathematical  models  of 
equipments  to  be  calculated.     PACER  as  supplied  by  Prof.  Shannon  con- 
tained forty- two  equipment  subroutines  in  the  PACER  equipment  subrou- 
tine library.     The  user  may  supply  such  additional  equipment  subroutines 
as  he  may  require. 

4.  Data  deck  -  the  data  deck  describes  the  system  or  process 
configuration  and  any  known  stream  variables  or  equipment  parameters. 
It  also  includes  certain  dimension  and  control  parameters  as  well  as 
certain  other  operational  indicators  or  controllers.    The  principal  con- 
stituents of  the  data  deck  are: 

(a)  Process  Matrix  -  the  process  matrix  is  the  principal 
segment  of  the  data  deck.     It  describes  the  sytem  or  problem  configuration 
and  indicates  stream  interrelationships.     Column  1  contains  the  equip- 
ment subroutine  which  performs  the  mathematical  calculation  of  the  equip- 
ment listed  in  column  1  „     The  columns  after  column  2  indicate  the  numbers 
of  the  streams  flowing  into  and  out  of  equipment  number  listed  in  column 
1.    All  input  streams  are  listed  first.    All  output  streams  must  be  preceded 
by  a  minus  sign. 

(b)  Equipment  Parameters  Matrix  -  contains  for  each  equipment 
listed  in  the  process  matrix  a  list  of  equipment  parameters  such  as  weight, 
volume,  power,  or  heat  transfer  coefficient.     Column  1  indicates  the 
equipment  number.    All  columns  after  column  1  contain  equipment  para- 
meters for  that  equipment.    The  equipment  parameters  are  characteristic 

of  that  equipment  only. 

(c)  Stream  Variables  Matrix  -  contains  stream  variables  for 
each  stream  listed  in  the  process  matrix.    The  stream  variables  represent 
such  characteristic  items  as  pressure,  temperature,  specific  heat,  total 
flow  rate,  component  flow  rates,  etc.     Stream  variables  are  consistent 
for  each  stream.     Column  1  contains  the  stream  number.     Column  2  is  a 
flag  indicating  the  status  of  the  stream.     If  the  flag  equals  1,   the  stream 


is  a  feed  stream.     If  it  equals  2,   it  is  a  product  stream.     If  it  is  equal 
to  0,  it  is  an  internal  stream.    The  remaining  columns  contain  the  values 
of  the  stream  variables  for  that  stream. 

(d)  Convergence  Criteria  -   contains  in  each  column  an  abso- 
lute fractional  convergence  value  for  each  corresponding  column  in  the 
stream  matrix.     Convergence  of  an  item  in  the  stream  variables  matrix 
occurs  when  the  value  calculated  for  that  variable  does  not  vary  from  the 
current  value  of  that  variable  stored  in  the  stream  variables  matrix  by 
more  than  the  absolute  convergence  criterion.    The  convergence  criterion 
for  a  particular  variable  location  (column)  is  the  same  for  all  streams. 

(e)  Calculation  Order  Listing  -  if  provided  or  specified  by  the 
user  it  indicates  the  order  in  which  the  equipment  will  be  calculated  as 
well  as  the  mode  (direct  or  iterative)  by  which  the  equipments  will  be 
calculated.     Column  1  indicates  the  equipment  number  and  column  2 

the  mode  of  calculation.    A  value  of  1  in  the  calculation  mode  column 
indicates  that  the  equipment  is  to  be  calculated  directly.    A  value  of  2 
indicates  that  the  equipment  is  to  be  calculated  iteratively.     If  the  cal- 
culation mode  is  3  the  equipment  is  the  last  equipment  in  a  list  of 
equipment  to  be  calculated  in  an  iterative  mode. 

(f)  Equipment  Subroutine  Listing  -  a  listing  which  establishes 
the  correspondence  between  the  subroutine  name  and  subroutine  number 
to  facilitate  calling  of  proper  subroutine  for  equipment  calculations. 

(g)  Other  control  and  informational  parameters  and  data  such 
as  problem  dimension  parameters,  problem  control  parameters ,  and  title 
of  simulation  run.    These  items  are  specifically  defined  in  Appendix  A. 
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Chapter  IV 

MODIFICATIONS  MADE  TO  PACER 

Modifications  Mode  to  PACER.    The  PACER  source  deck  supplied 
by  Prof.  Shannon  was  impressed  on  a  2400  foot  magnetic  iape  at  556 
bits/inch,   80  characters/record,  blocked  1.    With  the  assistance  of 
the  University  of  Kansas  Computation  Center,  the  magnetic  tape  was 
converted  into  an  interpreted,  punched  card  deck  of  approximately 
11,500  cards. 

Detailed  analysis  of  the  interpreted  deck  revealed  that: 

(1)  The  program  was  organized  into  eight  16K  links  as  well  as 
a  starter  link  and  several  sample  problems.    The  eight  links  were 
arranged  as  follows: 

a.  LINK  1  -  contained  mainline  segment  for  conducting  iterative 
calculations,  calling  and  testing  subroutines  and  certain  equipment 
subroutines  used  in  the  contact  sulfuric  acid  plant  simulation. 

b.  LINK  3  -  consisted  of  mainline  segment  which  directs  search 
for  list  of  equipment  when  order  of  calculation  is  not  specified,  sub- 
routines for  formulation  of  equipment  calculation  list,  and  a  subroutine 
for  recording  values  calculated  in  direct  calculations. 

c.  LINK  4  -  composed  of  a  mainline  segment  which  directed 
input/output  operations  and  subroutines  which  effected  input/output 
operations  and  stream  and  equipment  flagging. 

d.  LINKS  2,5,6,7,8  -  contained  remaining  equipment  and 
information  modification  and  gathering  subroutines  used  in  contact 
sulfuric  acid  plant. 

(2)  The  linking  procedure  employed  in  the  program  was  not 
compatible  with  the  GE  625  system  installed  at  the  University  of 
Kansas . 

(3)  The  program  was  written  in  FORTRAN  IIS 

Further  analysis  indicated  that  the  11,500  cards  could  be 
reduced  to  approximately  2  800  by  removing  all  those  cards  pertaining 
to  equipment  subroutines  used  in  the  contact  sulfuric-acid- plant 
simulation  as  well  as  those  used  in  the  sample  problems  at  the  end  of 
the  source  deck.    The  cards  remaining  after  this  modification  constituted 
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the  PACER  executive  routine.     However,  it  still  existed  in  a  non- 
compatible  linked  version  and  was  written  in  FORTRAN  II  which  was 
not  permitted  on  the  GE  625  system. 

Accordingly,  two  additional  modifications  were  required  to 
make  the  PACER  executive  routine  compatible  with  the  GE  625  system. 
First,  the  program  had  to  be  converted  to  a  linked  version  acceptable 
by  the  GE  62  5  and  the  program  had  to  be  rewritten  in  FORTRAN  IV.    A 
detailed  investigation  of  the  organization  of  the  executive  routine  and 
the  simulation  model  planned  for  development  indicated  that  the  PACER 
executive  and  the  simulation  model  could  be  combined  into  a  non- linked 
program  which  would  not  exceed  the  memory  available  on  the  625  . 

The  conversion  from  a  linked  to  a  non- linked  version  created 
several  exacting  problems.     First,  the  combination  of  the  linked  main- 
line segments  into  one  non- linked  mainline  program  created  a  statement 
number  redundancy.     Second,  it  caused  the  invalidation  of  a  major 
portion  of  the  linked  program  transfer  logic.    This  problem  required 
extensive  effort  to  resolve. 

The  rewriting  of  the  executive  routine  in  FORTRAN  W  required 
considerably  less  effort  than  the  conversion  to  a  non-linked  program. 
The  major  changes  required  were  the  rewriting  of  arithmetic  function 
statements  and  input/output  instructions. 
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Chapter  V 

DEVELOPMENT  OF  STEAM  SYSTEM  MODEL  AND 
TESTING  OF  THE  MODIFIED  PACER  EXECUTIVE  ROUTINE 

Objectives  of  Testing.    The  principal  objectives  of  testing 
conducted  with  the  modified  PACER  executive  routine  were: 

(1)  To  verify  the  validity  of  the  modifications  made  to  the 
original  version  of  PACER; 

(2)  To  demonstrate  the  capability  of  PACER  to  simply  and 
accurately  simulate  the  performance  of  varying  configurations  of  com- 
plex chemical  engineering  systems.     In  this  regard  it  should  be 
pointed  out  that  specific  care  was  taken  throughout  the  course  of  the 
modification  of  PACER  to  insure  that  no  alterations  were  made  in 
either  PACER  logic  or  method  of  solution.     Finally, 

(3)  To  provide  the  experimental  and  factual  basis  for  evaluating 
the  potential  usefulness  of  the  modified  PACER  executive  routine  in  a 
chemical  engineering  curriculum.    This  objective  was  facilitated  by 
the  steam  power  plant  model  developed  by  the  author  for  testing  with 
the  modified  PACER  program. 

Development  of  Simulation  Model.    Feasibility  or  validity  testing 
of  the  modified  PACER  program  presented  the  contrasting  alternatives 
of  (1)  simulation  of  a  segment  of  the  contact  sulfuric  acid  plant  by 
use  of  the  equipment  subroutines  supplied  by  Prof.  Shannon  or  (2) 
development  of  a  new  simulation  model.    The  second  alternative  was 
selected  since  it  provided  the  opportunity  for  development  of  a  more 
thorough  understanding  of  PACER  program  logic. 

The  simulation  model  selected  was  that  of  a  conventional 
steam  power  system  utilizing  fuel  oil  as  its  energy  source.    The  reasons 
supporting  this  selection  were: 

(1)     The  familiarity  of  the  author  with  steam  power  systems. 
This  familiarity  was  gained  by  the  author  from  service  for  two  years 
as  an  instructor  in  Thermodynamics  at  the  United  States  Naval 
Academy,  as  well  as  practical  acquaintance  with  steam  power  systems 
as  the  primary  means  for  naval  ship  propulsion. 


13 

(2)  The  general  applicability  of  steam  systems  to  chemical 
engineering . 

(3)  The  ready  adaptability  of  such  a  model  to  the  undergraduate 
chemical  engineering  curriculum  at  the  University  of  Kansas. 

Accordingly,  the  author  undertook  the  development  of  the 
essential  equipment  subroutines  necessary  for  the  approximation  of 
a  steam  power  system.    As  envisioned,  the  steam  power  system  would 
consist  of: 

(1)  A  Turbine 

(2)  A  Condenser 

(3)  A  Cooling  Water  Circulating  Pump 

(4)  A  Condensate  Pump 

(5)  A  Heater  Mixer 

(6)  A  Feed  Pump 

(7)  An  Economizer 

(8)  A  Boiler 

(9)  A  Flow  Splitter 

(10)  A  Superheater 

(11)  A  Fuel  Oil  Heater 

(12)  A  Forced  Draft  Blower 

(13)  A  Flow  Mixer 

Such  a  system,  it  was  envisioned,  would  be  amenable  to  demon- 
strating the  ability  of  PACER  to  perform  iterative  material  and  energy 
balances  by  the  multiple  application  of  individual  equipment  sub- 
routines. 

The  basic  system  as  envisioned  and  developed  by  the  author  is 
presented  in  Figure  3.    Figures  4  and  5  present  alternate  configurations 
of  the  basic  system.    Table  1  provides  a  stream  legend  for  Figures  3-5. 
Figure  6  presents  the  process  diagrams  for  the  basic  system  on 
enthalpy-- entropy  coordinates.    Analysis  and  development  of  mathe- 
matical models  for  individual  components  of  the  basic  system  are 
presented  in  Figures  7-19. 
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Stream  Number  Description 

1  Steam 

2  Exhaust  Steam 

3  Impurities 

4  Losses 

5  Condensate 

6  Cooling  Water 

7  Cooling  Water 

8  Cooling  Water 

9  Auxiliary  Steam 

10  .        Auxiliary  Steam 

11  Condensate 

12  Auxiliary  Steam 

13  Auxiliary  Steam 

14  Feed 

15  Feed 

16  Auxiliary  Steam 

17  Auxiliary  Steam 

18  Feed 

19  Steam 

20  Fuel  Oil 

21  Boiler  Fuel  Oil 

22  Boiler  Fuel  Oil 

23  Boiler  Blowdown  of  Impurities  and 

Accompanying  Feed 

24  Fuel  Oil 

25  Superheater  Fuel  Oil 

26  Superheater  Fuel  Oil 

27  Auxiliary  Steam 

28  Auxiliary  Steam 

29  Air 

30  Air 

31  Auxiliary  Steam 

32  Auxiliary  Steam 

Table  1.     Steam  Power  Plant  Stream  Legend 
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Stream  Number  Description 

33  Flue  Gases 

34  Steam 

35  Flue  Gases 

36  Make-Up  Feed 


Table  1 .     (continued) 
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Power  =  W  .  x  (h,   -  h0)  x  N 

steam  1         2  m 

h0  =  (h,   -  h,    )  x  N    . 
2       v   1         Is  ei 


W  .  -  steam  flow  rate 

steam 

h,  -  entering  enthalpy 

N    .  -  internal  engine  efficiency 

N  -  mechanical  efficiency 
m 

h,  -  isentropic  expansion  enthalpy 

ru  -  exhaust  enthalpy 


POWER 


Figure  7.    Turbine 
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Q 

steam        .'  ~i         5'  inpur 

W  =  Q/C        x    (T7  -  Tc) 

cw  cw        x   7         6 

W.  =  F,    x  W   . 

inpur  1  steam 

W,  =  F0   x  W   . 

loss  2  steam 

cond  steam  inpur  muf  loss 

THETAM  -  THETA1    -    THETA2 
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AREA  -  Q/(U  x  THETAM) 
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loss 

THETAM 
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W  . 
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cw 


rl 
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cona 


-  heat  transferred 

-  inpurities  flow  rate 

-  make-up  feed  flow  rate 

-  temperature 

-  losses  flow  factor 

-  water  flow  rate  in  blowdown 

-  system  losses 

-  log  mean  temperature  difference 

-  larger  temperature  difference 

-  smaller  temperature  difference 

-  overall  heat  transfer  area 

-  condenser  heat  transfer  area 

-  steam  flow  rate 

,h^fi   -  enthalpy 

-  cooling  water  specific  heat 

-  inpurities  flow  factor 

-  condensate  flow  rate 


Figure  8  .     Condenser 
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-  saturated  steam  flow  rate 

-  pumped  liquid  flow  rate 

-  specific  volume  of  pumped  fluid 

-  pressure 

h7,h,  n      -  enthalpy 

-  isentropic  expansion  enthalpy 

-  internal  engine  efficiency 

-  mechanical  efficiency 
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Figure  9.     Circulating  Pump 
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Figure  10.     Condensate  Pump 
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Figure  11.     Heater- Mixer 
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-  isentropic  expansion  enthalpy 
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Figure  12.     Feed  Pump 
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Figure  13.     Economizer 
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Figure  14.     Boiler 
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No  mathemaical  calculations  take  place  in 
the  splitter.  The  splitter  merely  serves  as 
a  convenient  manner  for  splitting  the  sat- 
urated steam  flov/  produced  in  the  boiler  into 
its  various  components  v/ithout  needlessly  en- 
cumbering the  boiler  calculations. 


Figure  15.  Splitter 
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Figure  16.     Superheater 
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Figure  17.    Fuel  Oil  Heater 
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Figure  18.     Forced  Draft  Blower 
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Plan  of  Testing.    The  development  of  a  series  of  approximations  to 
achieve  the  objectives  established  above  required  a  systematic  analysis 
to  determine  if  any  of  the  objectives  were  complementary  or  interdepend- 
ent, and,  if  so,  in  what  manner  and  in  what  order. 

An  analysis  of  these  relationships  revealed  that  the  objective  of 
verifying  the  validity  of  the  modifications  made  to  PACER  was  the  objec- 
tive on  which  all  the  other  objectives  depended.    Additionally,  it  was 
concluded  that  objectives  two  and  three  could  be  achieved  as  derivatives 
of  objective  one.    That  is,  if  simulation  tests  could  be  developed  which 
would  establish  the  validity  of  the  modifications  and  the  accuracy  of  the 
model,  then  the  efficacy  of  PACER  and  the  model  and  their  usefulness  in 
chemical  engineering  education  could  be  determined  by  derivation  from 
those  results . 

Accordingly,  a  series  of  six  simulations  were  scheduled  for  testing 
with  the  modified  PACER  program.    The  first,   Run  01,  was  designed  to 
approximate  the  performance  of  the  basic  system  with  the  order  of  calcu- 
lations specified  by  the  user.    The  second  simulation",  designated  Run  03, 
involved  the  same  basic  system  but  did  not  specify  the  order  of  calcula- 
tions.    The  accuracy  of  the  solutions  was  not  of  primary  concern  in  these 
tests.     Rather,  the  objective  was  to  ascertain  if  the  modified  program  were 
processing  the  data  in  accordance  with  the  methods  described  in  the  docu- 
mentation provided  by  Prof,  Shannon  and  inferred  from  an  analysis  of  the 
original  PACER. 

Two  additional  system  configurations  were  tested  to  demonstrate  the 
ability  of  PACER  to  provide  accurate  solutions  of  varying  configurations  of 
the  same  system.    The  first  of  the  altered  configurations,  labeled  Run  02, 
involved  all  the  components  of  the  basic  system  except  the  superheater. 
The  second  altered  configuration  consisted  of  the  basic  system  less  the 
economizer.    This  run  was  designated  Run  04.    Both  Runs  02  and  04  were 
run  with  the  order  of  calculation  specified  and  unspecified.    The  results 
obtained  from  these  runs  were  compared  with  hand- calculated  values  and 
corrections  as  necessary,  and  were  made  in  the  models  until  the  computed 
values  agreed  with  the  hand  calcuations  for  all  cases. 

The  final  results  of  all  runs  are  presented  in  Appendix  B,  with 
extensive  diagnostic  print-out  (a  PACER  option)  provided  for  Runs  01  and 
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03.    A  guide  for  the  interpretation  of  the  results  of  the  simulations  accom- 
panies the  results  in  Appendix  B. 
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Chapter  VI 

ANALYSIS  OF  TEST  RESULTS 

Accuracy  of  Results.    Acceptable  limits  of  accuracy  vary  with  the 
nature  of  the  problem  being  evaluated.    Quite  clearly,  the  acceptable 
limits  of  accuracy  in  determining  the  configuration  of  the  critical  mass 
of  a  fissionable  material  would  be  much  more  rigorous  than  those  to  be 
imposed  on  the  sytem  parameters  and  configuration  of  a  steam  power 
plant.     In  the  latter  case,   slide  rule  accuracy  is  normally  acceptable. 
Accordingly,  a  convergence  criterion  of  .02  (2%) — an  approximation  of 
slide  rule  accuracy — was  selected  for  all  steam  variables  other  than  those 
non- varying  items  such  as  stream  number  and  flag  value.    With  this  cri- 
terion, all  runs  converged  to  solution  in  no  more  than  five  iterations. 

The  values  of  stream  variable  calculated  by  PACER  were  compared 
to  hand-calculated  values  and  in  all  cases  were  within  slide  rule  accuracy 
of  hand-calculated  values. 

Effect  of  Specifying  Order  of  Calculations.     PACER  undertakes  trial 
and  error  (iterative)  solutions  in  all  cases  where  the  order  of  calculations 
is  not  specified  or  where  an  iterative  calculation  is  specified.     Under 
these  conditons,  the  advantage  of  specifying  the  order  of  equipment  cal- 
culations was  not  clear.     In  Run  02   (without  superheating)  both  the  calcu- 
lations with  and  without  the  order  of  calculations  specified  converged  in 
five  loops.    In  Run  04  (without  economizer)  the  calculations  with  the  order 
of  calculations  specified  converged  in  three  iterative  loops,  whereas  the 
calculations  without  the  order  of  calculations  specified  converged  in  four 
loops.     Runs  01  and  03  (basic  system)  converged  in  five  and  four  loops 
respectively. 

The  inconclusive  nature  of  these  results  would  seem  to  indicate 
that  the  PACER  method  of  determination  of  order  of  calculation  is  as  effec- 
tive in  promoting  convergence  as  intuitive  designation  of  the  order  of 
calculation — or  at  least  no  more  ineffective. 

Another  aspect  of  this  consideration  is  the  comparison  of  the  values 
of  stream  variables  obtained  with  different  orders  of  calculations.    Quite 
clearly  convergence  values  will  be  expected  to  vary  depending  on  the 
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logical  location  of  equipment  in  the  order  of  calculations.    While  this 
expectation  was  confirmed,  differences  were  in  all  cases  less  than  the 
convergence  criterion. 

Effects  of  Changes  in  System  Configuration.  Two  major  changes 
to  the  basic  steam  power  system  were  simulated  with  the  modified  PACER 
program.  The  first  change  involved  removal  of  the  superheater  from  the 
basic  system.  This  change  necessitated  the  use  of  saturated  steam  in 
the  turbine  with  the  result  that  for  the  same  turbine  power  requirement, 
steam  flow  increased  12%.  However,  fuel  oil  flow  increased  only  1%. 
Condenser  surface  area  increased  6%. 

The  second  change  involved  removal  of  the  economizer  from  the 
steam  power  system.    This  change  necessitated  the  addition  of   the  heat 
that  would  have  been  added  by  flue  gases  in  the  economizer  as  the  heat 
of  combustion  of  additional  fuel  oil.    This  change  resulted  in  a  20% 
increase  in  fuel  oil  flow. 

Realism  of  the  Models .    The  qualitative  and  quantitative  nature  of 
the  stream  variables  tended  in  all  cases  to  verify  the  realistic  nature  of 
the  models  developed  by  the  author.    For  example,  the  use  of  saturated 
steam  in  lieu  of  superheated  steam  would  undoubtedly  require  more  steam 
to  produce  the  same  turbine  power.    The  greater  steam  flow  would  norm- 
ally require  more  fuel  oil.     However,  the  fact  that  the  additional  steam 
was  saturated  and  the  original  superheating  load  was  no  longer  required 
had  an  offsetting  effect  on  fuel  oil  flow. 

Additionally,  as  was  demonstrated  in  Run  04,  the  loss  of  the  regen- 
erative heating  affect  of  the  economizer  would  quite  clearly  require  more 
fuel  oil  for  the  same  power  requirement. 

While  these  results  tend  to  verify  the  realistic  nature  of  the  steam 
power  system  developed  by  the  author,  it  should  be  understood  that  no 
attempt  was  made  in  the  development  of  the  model  to  completely  depict 
reality.  Rather  an  attempt  was  made  to  provide  a  semi-realistic  model 
v/hich  could  be  refined  subsequently.  It  should  be  noted  in  this  regard 
that  most  of  the  data  and  mathematical  bases  for  equipment  and  sytems 
models  are  based  on  reliable  data  developed  from  many  years  in  the 
design,  development  and  operation  of  naval  steam  power  plants. 
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Validity  of  Modifications  Made  to  PACER.    The  close  conformance 
of  the  values  obtained  in  the  simulations  with  hand-calculated  values, 
as  well  as  the  general  and  specific  correspondence  of  results  to  format 
and  arrangement  of  data  as  described  in  documentation  provided  by  Prof. 
Shannon  clearly  establishes  the  validity  of  modifications  to  the  original 
PACER  program  o 


Chapter  VII 
AN  ASSESSMENT  OF  PACER 

General.    The  author  has  demonstrated  the  concept,  principles  and 
method  of  application  of  the  PACER  simulation  program  by  the  develop- 
ment of  a  complex  model  of  a  steam  power  system,  and  successfully 
approximating  its  performance  on  a  modified  version  of  the  original  PACER. 
The  validity  of  the  modifications  made  to  PACER  have  been  conclusively 
established  by  the  accuracy  of  the  results  achieved.    The  realistic  nature 
of  the  model  has  been  tentatively  established  by  the  qualitative  and 
quantitative  nature  of  the  results. 

PACER  As  A  Simulation  Device.     PACER  has  been  demonstrated  to  be 
a  powerful,  flexible  and  efficient   method  for  simulation  of  complex  con- 
figurations of  chemical  engineering  systems. 

However,  certain  areas  where  improvements  in  PACER  method  of 
operations  or  concept  could  be  made  were  noted  and  are  discussed  here 
for  edification  of  potential  users. 

(1)  The  addition  or  removal  of  equipment  subroutines  from  PACER 
is  cumbersome  with  present  arrangement  of  calling  segments  and  sub- 
routines and  LPHA  listing.    A  more  permanent,  non-mnenonic  name  for 
equipment  subroutines  requiring  only  a  change  in  subroutine  name  to  con- 
form to  permanently  established  but  unused  name  is  one  method  of  allevi- 
ating this  problem.    This  would  obviate  changes  in  mainline  program  which 
become  more  difficult  when  executive  program  is  on  tape. 

(2)  Consideration  should  be  given  to  establishment  of  permanent 
STRMI  and  STRMO  matrices  to  accommodate  all  streams  on  a  permanent 
basis  for  a  run.    This  would  facilitate  multiple  calling  of  stream  variable 
and  information  transmittal  in  a  particular  run.     However,  it  could  develop 
that  the  ultimate  result  would  be  to  trade  the  present  inflexibility  of  a 
temporary  STRMI  and  STRMO  matrix  for  an  inflexible  permanent  matrix. 

(3)  Consideration  should  be  made  for  incorporating  capital  invest- 
ment costs  into  equipment  subroutines .    This  information  could  take  the 
form  of  a  cost  which  varies  with  heat  transfer  area,  pressure,  temperature, 
etc.     Incorporation  of  this  data  into  the  basic  structure  of  an  equipment 


subroutine  is  a  relatively  simple  matter — primarily  a  function  of  availa- 
bility of  cost  information  and  its  relationship  to  design  or  operating 
parameters . 

Use  of  PACER  In  A  Chemical  Engineering  Curriculum.     In  evaluating 
the  results  of  the  simulation  of  the  C-I-L  contact  sulfuric  acid  plant 
using  the  PACER  program /   Prof.  Shannon  remarked: 

It  should  be  emphasized  that  this  study,  although  it 
involved  a  unique  and  very  close  cooperation  with 
industry,  was  in  the  main  an  exercise  in  the  education 
of  engineering  students.    The  students  were  concerned 
with  the  detailed  modelling  of  specific  equipment,  yet 
had  a  very  broad  exposure  to  a  real  industrial  plant  and 
its  over- all  behavior.    They  learned  of  the  difficulties 
and  importance  of  obtaining  good  plant  data,  and  of 
writing  successful  programs  of  greater  complexity  than 
previously  encountered.    The  project  was  admittedly  a 
very  ambitious  one;  however  the  large  effort  expended 
has  placed  PACER  into  the  hands  of  the  staff  for  future 
use  in  the  undergraduate  and  graduate  program. 

The  staff  generally  agree  that  most  of  the  students  have 
benefitted  by  this  exposure  to  the  quite  new  concept  of        ?7 
information  handling  and  decision  making  utilizing  PACER. 

PACER  then  has  already  demonstrated  its  effectiveness  in  instilling 
in  students  some  aspects  cf  "engineering  experience"  which  would  not 
normally  be  obtainable  outside  an  industrial  environment.    This,  as 
stated  in  the  University  of  Michigan  Report  on  the  uses  of  computers  in 
engineering  education,  is  the  raison  d'etre  of  engineering  design  courses. 

It  is  further  considered  that  PACER  can  be  an  effective  means  of 
developing  a  systematic  approach  to  problem  solving  by  chemical  engi- 
neering students.    The  demands  of  computer  programming  in  development 
of  equipment  subroutines  and  systems  concepts  requires  strict  attention 
to  problem  structuring  and  formulation. 

Recommmendations .     In  consideration  of  the  results  of  the  simu- 
lations conducted  with  the  modified  PACER  program  and  the  conclusions 
drawn  above,  it  is  recommended  that: 

(1)     Consideration  be  given  to  the  incorporation  of  the  modified 
PACER  program  into  existing  undergraduate  courses  involving  basic  prob- 
lems in  chemical  engineering.    The  steam   power  system  as  developed  by 
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the  author  is  amenable  to  immediate  use  on  an  individual  or  term  project 

basis/  and 

(2)  After  suitable  evaluation  of  the  effectiveness  of  the  initial  use 
of  PACER  that  consideration  be  given  to  further  and  more  extensive  use  of 
PACER  in  graduate  and  research  projects  at  the  University  of  Kansas. 
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Appendix  A 

USER'S  GUIDE  FOR  THE  MODIFIED  PACER  SIMULATION  PROGRAM 


General,    The  most  significant  feature  of  the  PACER  simulation 
method  is  its  ability  to  perform  realistic  and  accurate  simulations  of 
varying  configurations  of  complex  chemical  engineering  systems  with  a 
minimal  effort  on  the  part  of  the  user.    The  element  of  PACER  which 
makes  this  capability  possible  is  a  highly- structured  data  deck  which 
permits  the  user  to  fully  characterize  system  configuration  and  the  re- 
straints or  parameters  under  which  it  operates. 

Successful  use  of  PACER  is' dependent  upon  strict  adherence  to  a 
prescribed  set  of  instructions  for  the  formulation,  assembly  and  modifi- 
cation of  the  data  deck.     It  is  the  purpose  of  this  appendix  to  provide 
the  necessary  instructions  in  a  manner  which  will  permit  a  user  to 
employ  PACER' s  powerful  calculational  skills  without  having  to  master 
either  the  logic  or  structure  of  PACER.     One  exception  exists  to  this 
condition.    That  occurs  when  it  is  necessary  to  develop  or  modify  equip- 
ment subroutines  or  add  them  to  or  delete  them  from  the  program.     In  that 
case  a  minimal  knowledge  of  the  structure  and  logic  of  PACER  is  required 
to  insure  the  subroutines  conform  to  PACER  logic  and  that  they  are  properly 
entered  into  PACER  structure.     Both  of  these  aspects  will  be  treated  in 
this  appendix. 

System  Conceptualization.    The  initial  step  in  preparing  a  system 
configuration  for  simulation  by  PACER  is  to  establish  the  system  concept 
and  commit  the  concept  to  formalization  by  preparation  of  a  flow  diagram. 
This  flow  diagram  should  break  equipment  assemblies  into  subassemblies 
and  equipments  capable  of  being  mathematically  modelled.     For  example, 
a  boiler  normally  contains  an  economizer  and  a  superheater  as  integral 
equipments,  and  when  one  normally  refers  to  a  boiler  it  is  assumed  that 
he  includes  the  economizer  and  superheater  if  they  are  provided.     However, 
breaking  these  elements  out  and  handling  them  as  separate  equipments 
facilitates  mathematical  and  system  modelling  as  well  as  increasing  the 
flexibility  of  the  model.     Further,  the  functional  components  if  of  suffi- 
cient general  nature  could  be  used  in  multiple  applications  within  a 
single  simulation  or  as  a  related  equipment  in  subsequent  simulations. 
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In  developing  the  system  flow  diagram  all  equipments ,  the  per- 
formance of  which  will  be  simulated  by  mathematical  models  (equipment 
subroutines),   should  be  numbered.    No  redundancy  in  equipment  number- 
ing is  permitted  although  equipment  numbers  need  not  be  successive. 
While  equipments  are  being  numbered,  they  should  also  by  named.    For 
convenience  it  would  be  preferable  to  assign  them  a  six-character  alpha- 
numeric name  which  is  mnemonically  related  to  the  equipment  function, 
e.g.  TURBN5  for  turbine  and  CNDNR5  for  condenser.    No  redundancy  is 
permitted  in  the  last  3  characters  of  a  unit  name,  although  individual 
characters  within  a  unit  name  may  be  repeated.    For  example,  TURBN5 
and  TURBN6  are  acceptable,  while  SPLTR5  and  SPHTR5    are  not. 

After  all  equipments  have  been  numbered  and  named,  all  inter- 
connecting streams  should  be  numbered.     No  redundancy  is  permitted 
in  stream  numbers.     However,   stream  numbers  need  not  be  successive 
nor  is  redundancy  between  stream  and  equipment  numbers  prohibited. 

Equipment  Subroutine  Format.    The  equipment  subroutine  is  the  one 
facet  of  PACER  that  poses  a  programming  burden  on  the  user.    If  an 
equipment  subroutine  which  satisfies  the  analytical  and  calculational 
requirements  of  an  equipment  which  is  to  be  simulated  is  not  available 
in  the  PACER  subroutine  library,  the  user  must  develop  an  equipment 
subroutine  which  accomplishes  these  objectives. 

Several  challenges  confront  the  user  developing  an  equipment  sub- 
routine for  use  in  a  PACER  simulation.    The  first  challenge  is  to  insure 
that  the  equipment  to  be  modelled  is,  in  fact,  the  equipment  which  should 
be  modelled.     In  this  respect  a  general  effort  should  be  made  to  break 
equipment  assemblies  down  into  their  functional  elements,  that  is,  as 
previously  mentioned,  a  boiler  into  a  boiler,  an  economizer  and  a  super- 
heater.    Mathematical  analysis  is  facilitated  and  simplified  by  this 
approach.     Further,   specific  effort  should  be  directed  toward  incorporating 
generality  into  subroutines.    That  is,  an  economizer  is  in  essence  a  two- 
fluid  heat  exchanger.     If  the  equipment  subroutine  can  be  written  for  a 
general  two- fluid  heat  exchanger  which  is  in  this  application  called  an 
economizer,  then  the  equipment  subroutine  can  be  used  in  subsequent 
future  applications  where  a  two- fluid  heat  exchanger  is  employed.    The 
fact  of  this  generality  should  be  documented  into  the  subroutine  by 
appropriate  comments,, 
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Once  the  mathematical  representation  of  an  equipment  has  been 
developed  it  is  necessary  to  translate  this  into  the  syntax  employed  by 
PACER.    This  task  can  be  mastered  quickly  by  knowledge  of  a  few  basic 
concepts  of  PACER. 

PACER  performs  equipment  calculations  by  using  mathematical 
relationships  provided  by  the  user  in  the  equipment   subroutine  and 
stream  and  equipment  data  provided  by  the  user  in  the  stream  (SN)  and 
equipment  (EN)  matrices  of  the  data  deck  (to  be  discussed  further  in 
section  dealing  with  the  data  deck).     PACER  calls  data  out  of  the  stream 
matrix  into  temporary  STRMI  (inlet  streams)  and  STRMO  (outlet  streams) 
for  the  equipment  to  be  calculated.    The  inlet  and  outlet  streams  for  a 
particular  equipment  are  designated  by  the  user  in  the  process  matrix 
(KPM) ,     Stream  variable  information  is  then  ready  for  use  in  PACER  cal- 
culations.    Equipment  parameters  are  called  from  and  returned  directly 
to  the  EN  matrix  for  the  equipment  to  be  calculated.    Once  the  equipment 
has  been  calculated  the  values  in  the  STRMI  and  STRMO  matrices  are 
returned  to  permanent  storage  in  the  SN  matrix.    The  STRMI  and  STRMO 
matrices  are  then  available  for  temporary  storage  of  inlet  and  outlet 
stream  variables  for  next  equipment  to  be  calculated. 

It  is  then  incumbent  upon  the  user  to  transcribe  his  mathematical 
notation  into  a  notation  which  indicates  to  PACER  where  in  the  STRMI/ 
STRMO  or  EN  matrices  the  variable  or  parameter  to  be  used  is  to  be  lo- 
cated.   For  example ,  let  us  convert  the  mathematical  equation  for  deter- 
mining turbine  stream  flow  rate  into  a  notation  usable  by  PACER.    The 
equation  (neglecting  dimensional  constants)  is: 

W  =   P  4-Wk 

In  FORTRAN  this  is  represented  as 

W  =    P/WK 

However,  in  PACER  notation  it  would  be  written  as 

STRMI  (1,10)  =  EN  (NE,3)  /  (STRMI(1,5)  -  STRMO(l,5)) 
(where  WK  -  (STRMI(1,5)  -  STRMO(l,5))) 

In  the  above  notation 

STRMI(1,10)  -  indicates  the  tenth  variable  in  the  first  inlet  stream 
to  the  equipment  being  calculated.     It  also  represents  the  tenth 
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variable  in  the  stream  matrix  vector  for  the  stream  which  is  the 
first  input  stream  to  the  equipment.     In  the  context  of  this  sample, 
this  variable  represents  the  stream  flow  for  that  stream. 

EN  (NE,3)  -  indicates  the  third  parameter  for  equipment  NE.     Equip- 
ment NE  is  the  equipment  number  for  the  equipment  being  calculated 
and  will  be  set  elsewhere  in  PACER.    The  turbine  power  rating  for 
equipment  NE  will  be  stored  in  this  location. 

STRMI  (1,5)  and  STRMO  (1,5)  -  represent  the  fifth  variable  in  the 
first  inlet  and  outlet  streams  respectively.  In  this  example  these 
variables  take  on  the  value  of  the  entering  and  leaving  enthalpies 
of  the  first  inlet  and  outlet  streams  to  equipment  NE. 

The  effect  of  the  above  calculation  is  to  calculate  a  value  for 
steam  flow  for  the  stream  number  which  is  the  first  inlet  stream  to  equip- 
ment NE.    After  calculation  this  calculated  value  is  stored  in  the  SN 
matrix  releasing  the  STRMI  and  STRMO  matrices  available  for  further 
calculations.     Since  the  values  of  EN(NE,3)  and  STRMI(1,5)  and  STRMO 
(1,5)  were  not  modified,  the  corresponding  values  in  the  EN  and  SN 
matrices  will  not  be  modified. 

Figures  Al.  and  A2  .  are  the  equipment  subroutines  for  the  mathe- 
matical models  developed  for  the  turbine  and  condenser  in  Chapter  4  (See 
Figures  7  and  8).     It  should  be  pointed  out  that  the  assignment  of  variable 
and  parameter  locations  in  the  SN,   EN,   STRMI  and  STRMO  matrices  are 
the  responsibility  of  the  user.     Further,   stream  variable  location  must 
be  consistent  for  all  streams  in  the  SN,   STRMI  and  STRMO  matrices. 
That  is  location  5  in  the  SN,   STRMI  and  STRMO  matrices  will   throughout 
the  course  of  a  run  represent  the  same  variable,  but  not  necessarily  the 
same  value.     For  example,   in  the  illustrative  example  discussed  above, 
location  5  represented  the  enthalpy  for  all  streams  but  the  enthalpy  for 
all  streams  was  not  necessarily  the  same  value  and  some  streams  did  not 
even  store  a  value  for  enthalpy  if  it  were  not  relevant  to  the  calculation 
of  that  equipment.     In  such  a  case  PACER  normally  assigns  a  value  of 
zero  to  that  location. 

One  additional  note  of  caution  is  that  a  common  and  dimension 
deck  identical  to  that  used  in  the  mainline  program  must  be  the  first  item 
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Figure  Al.  Subroutine  TURBNS 
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Figure   A2.    Contd. 
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in  each  subroutine.    This  facilitates  communication  of  information  from 
subroutines  to  mainline  and  other  subroutines. 

Incorporation  or  Removal  of  Equipment  Subroutines  into  PACER 
Structure.    The  flexible  nature  of  PACER  permits  the  incorporation  or 
removal  of  equipment  subroutines  to  conform  to  core  memory  limitations, 
to  permit  the  simulation  of  completely  different  systems  or  to  permit 
modification  of  existing  systems.     Such  flexibility  is  not  achieved,  how- 
ever, without  a  commensurate  demand  on  the  user.     In  this  case,   it  is 
the  requirement  that  the  incorporation  or  removal  of  a  subroutine  requires 
a  change  in  the  semi-permanent  structure  of  PACER. 

Each  subroutine,   in  addition  to  being  assigned  a  six-character, 
alpha-numeric  title,  is  assigned  an  equipment  subroutine  number  depending 
on  the  location  of  its  calling  statement  in  PACER.     Equipment  subroutines 
1-10  are  called  from  EQCAL  and  subroutines  11-30  are  called  from  the 
mainline  program.    The  location  of  these  calling  statements  and  maximum 
number  of  subroutines  are  a  function  of  the  system  to  be  simulated.     How- 
ever, a  note  of  caution.    An  increase  or  decrease  in  the  maximum  number 
of  equipment  subroutines  or  any  other  dimension  establishing  parameter 
requires  a  change  in  the  appropriate  matrix  or  vector  in  the  common  and 
dimension  deck. 

The  incorporation  of  a  subroutine  requires 

(1)  Obtaining  an  equipment  subroutine  number  from  those  not 
currently  being  used. 

(2)  Writing  a  calling  statement  for  that  subroutine. 

(3)  Providing  appropriate  transfer  directions  or  transfer  indicators. 

For  example,  the  addition  of  an  equipment  subroutine  COOLR8  as 
subroutine  number  5  would  require  the  insertion  of  the  following  FORTRAN 
cards  in  EQCAL  immediately  following  statement  number  4  and  its  transfer 
card,  if  any,  and  immediately  preceding  statement  number  6: 

Column  No.   5         7 

5        CALL  COOLR8 
RETURN 

The  insertion  of  the  same  subroutine  as  equipment  subroutine  number  15 
would  require  the  insertion  of  the  following  cards  in  the  mainline  program 
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immediately  following  statement  number  14  and  its  following  transfer 
card: 

Column  No.        4       7 


15     CALL  COOLR8 
GO  TO  7777 

Or  the  addition  of  COOLR8  as  equipment  subroutine  number  21  or  any 
number  greater  than  21  would  require  the  insertion  of  the  following  cards 
in  the  proper  location  following  statement  number  52  01: 

Column  No,        4       7 


21     CALL  COO LR8 
GO  TO  7777 

The  removal  of  an  equipment  subroutine  requires  the  reverse  procedure — 
removal  of  the  calling  statement  and  the  transfer  card  and  replacement  of 
these  cards  with  the  transfer  card  assigned  the  statement  number  of  the 
deleted  subroutine.     For  example,  if  equipment  subroutine  15  were  deleted 
the  FORTRAN  cards  described  above  would  be  replaced  by  the  following 
FORTRAN  statement: 

Column  No.        4       7 


15     GO  TO  7777 

One  additional  change  is  required  in  addition  to  those  decribed 
above  when  an  equipment  subroutine  is  incorporated  into  or  removed  from 
PACERo    This  change  is  necessary  due  to  a  feature  of  PACER  which  estab- 
lishes the  correlation  between  equipment  subroutine  names  and  numbers 
and  permits  the  multiple  calling  of  equipment  subroutines  for  different 
equipments.    That  is  equipments  number  5  and  10  may  both  be  pumps 
which  can  be  calculated  by  subroutine  PUMP5.    This  feature  referred  to  as 
the  LPHA  listing  stores  the  last  three  characters  of  the  equipment  sub- 
routine name  in  the  position  corresponding  to  its  equipment  subroutine 
number.    The  format  of  the  LPHA  listing,  12/  (15A3,   10A3),  permits  the 
inclusion  of  the  first  25  equipments  on  card  i.     If  30  equipments  were 
used  two  cards  would  be  required--25  on  the  first  card  and  5  in  the  first 
15  spaces  of  the  second  card.    If  an  equipment  subroutine  number  is  not 
occupied,  its  three  characters  are  left  blank.    A  control  card  indicating 
the  number  of  equipments  for  which  an  LPHA  designation  should  be  read 
precedes  the  LPHA  listing  and  is  indicated  in  the  format  statement  as  an 
12  format. 
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If  equipment  subroutine  COOLR8  were  assigned  subroutine  number 
15,  the  following  entry  would  be  made  in  column  43-45  of  the  first  card 
of  the  LPHA  listing: 

Column  No.      43  -  45 


L  R  8 

If  equipment  subroutine  PUMPS  were  designated  equipment  number  28  the 
following  entry  would  be  made  in  columns  10-12  of  the  second  card  of  the 
LPHA  listing: 

Column  No.      10-12 


P  5 

If,  on  the  other  hand,  equipment  subroutine  PUMP5  were  removed  from  the 
PACER  structure,  columns  10-12  of  the  second  card  would  be  left  blank. 

Data  Deck.    The  data  deck  is  the  means  by  which  the  user  conveys 
to  the  PACER  executive  routine  is  a  detailed  discription  of  the  system  to 
be  simulated.     Only  one  simulation  run  may  be  processed  for  each  data 
deck  supplied.     However,   several  simulation  runs  may  be  processed 
consecutively  by  providing  consecutive  data  decks  for  each  run. 

In  general,  a  data  deck  for  a  simulation  run  is  composed  of  15 
groups  of  data.     However,  the  first  group,  the  LPHA  listing,   is  provided 
only  with  the  first  data  deck  of  a  multi-deck  execution.    All  decks  after 
the  first  deck  consist  of  only  14  data  groups.    The  data  groups  and  their 
required  order  of  assembly  are  listed   and  discussed  below.     In  addition, 
data  decks  for  the  six  simulations  run  to  test  the  feasibility  of  the  modified 
PACER  program  are  provided  as  illustrative  examples  at  the  end  of  this 
appendix.     (See  Figures  A3.  -  A8.). 

LPHA  Listing.     Format:    12/(1 5A3  ,  10A3) .    This  data  group  indicates 
the  allowable  names  of  PACER  equipment  subroutines  and  provides  the 
requisite  correspondence  between  the  equipment  subroutine  number,  the 
equipment  subroutine  name  and  the  equipment  number.    The  first  card 
contains  in  12  format  the  number  of  equipment  names  to  be  read  from  the 
LPHA  listing.    The  remaining     cards  contain  the  LPHA  listing  of  the  last 
three  characters  of  the  allowable  equipment  subroutines.    No  redundancy 
is  permitted  in  the  three- character  names  although  individual  character 
redundancy  is  permitted.    The  absolute  location  of  an  equipment  subroutine 
in  the  LPHA  listing  must  coincide  with  the  number  assigned  to  the  equipment 
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subroutine.    That  is,  the  equipment  subroutine  name  placed  in  the  25th 
location  in  the  LPHA  listing  must  be  the  last  three  characters  of  equipment 
number  25  . 

Title  Card.     Format:    2(12A3).    A  single  card  provided  for  the  assign- 
ment of  a  unique  identifying  title  to  each  run. 

Control  Parameters.  Format:  1015.  A  list  of  8  parameters  which 
control  the  processing  of  PACER  simulations.  The  parameters,  in  their 
order  of  listing  on  the  input  record,  are: 

1.  KRUN  -  Run  number. 

2.  KSETS  -  Controls  intermediate  printing  of  PACER  status. 
Following  is  a  listing  of  allowable  values  and  options  each 
provides . 

Value  Option 

0  Prints  initial  conditions  of  matrices,  inter- 

mediate conditions  of  matrices  every  KPRINT 
loops  and  final  results  of  run. 

2  In  addition  to  above,  prints  values  in  STRMI 
and  STRMO  after  each  calculation  and  number 
of  stream  variable  deviations  from  convergence 
criteria . 

3  In  addition  to  above,  prints  flag  lists  for 
streams  and  equipments. 

3.  KCLEAN  -  If  set  equal  to  zero,  it  causes  PACER  to  zero 
all  matrices  and  vectors  before  reading  in  data. 

4.  LOOPS  -  Maximum  number  of  loops  to  be  permitted  in 
iterative  calculations. 

5.  NOGO  -  Controls  simulation  options  if  convergence  is 
not  obtained  in  number  of  LOOPS  permitted .  If  NOGO  =  0, 
the  answers  obtained  during  the  last  loop  will  be  taken  as 
correct  and  the  calculation  will  continue.  If  NOGO  7^  0, 
additional  trial  and  error  testing  of  assumed  input  streams 
will  be  made  until  either  convergence  is  obtained,  execution 
time  expires  or  calculation  is  determined  to  be  impossible. 

6.  KARDS  -  If  greater  than  zero,  provides  a  deck  of  cards 
punched  with  the  stream  variables  matrix. 

7.  KPRINT  -  Controls  number  of  loops  between  printout  of 
stream  variables  matrix. 

8.  KONV  -  Not  used.     Leave  blank. 

Dimension  Parameters.     Format:    1015.    This  data  group  contains 

nine  parameters  which  control  the  allowable  dimensions  of  a  run.    Any 

change  in  these  variables  requires  an  identical  change  to  the  parameter 

where  it  appears  in  the  dimension  deck.    The  following  is  a  listing  of  the 

nine  parameters  in  the  order  of  their  recording  on  the  input  record: 

1.      N2MAX  -  maximum  number  of  rows  permitted  in  the  process 
matrix. 
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2.  N3MAX  -  maximum  number  of  columns  permitted  in  process 
matrix. 

3.  NEMAX  -  maximum  number  of  equipments  permitted  in 
system  to  be  simulated. 

4.  NSMAX  -  maximum  number  of  streams  permitted  in  system 
to  be  simulated. 

5 .  NELMAX  -  maximum  number  of  columns  permitted  in  EN 
matrix. 

6.  NECLMX  -  maximum  number  of  columns  permitted  in  equip- 
ment control  matrix. 

7.  NSLMAX  -  maximum  number  of  columns  permitted  in  SN 
matrix. 

8.  NSCLMX  -  maximum  number  of  columns  permitted  in  stream 
variables  control  matrix. 

9.  NOCOMP  -  maximum  number  of  components  in  each  stream. 
Not  used  other  than  for  information  purposes. 

Process  Matrix  (KPM) .    Format:    15,   1X2A3,   1X2A3,   1X915.    The 
process  matrix  contains  the  topology  or  configuration  of  the  system  to  be 
simulated.     It  is  preceded   by  a  single  card  (Format  15)  which  indicates 
the  number  of  rows  to  be  read  in  the  process  matrix.     Each  row  corresponds 
to  an  equipment  in  the  system  configuration.     Equipment  numbers  need  not 
be  consecutive  but  must  be  in  ascending  numerical  order.    The  first  item 
to  be  read  from  the  input  record  (Format  15)  is  the  equipment  number.    The 
second  and  third  items  are  the  first  three  and  second  three  characters, 
respectively,  of  the  6  character  equipment  subroutine  used  to  calculate 
that  equipment.    The  next  9  items  are  the  input  and  output  streams  for  that 
equipment.     Input  streams  are  listed  first  and  output  streams,  preceded  by 
a  negative  sign,  follow.    A  word  of  caution  is  necessary  at  this  point.    One 
must  remember  that  equipment  subroutines  are  written  in  anticipation  of  or 
in  conjunction  with  the  order  of  input  and  output  streams.    Accordingly, 
any  changes  in  the  order  of  input  and  output  streams  for  an  equipment  must 
be  checked  against  the  equipment  subroutine  to  insure  that  revised  input/ 
output  stream  listing  is  consistent  with  equipment  subroutine.     If  a  dis- 
parity exists  either  the  equipment  subroutine  or  the  order  of  input/output 
streams  must  be  revised  to  bring  them  into  conformity  with  each  other. 

Equipment  Parameter  Matrix  (EN).     Format:    5F15.5.    The  EN  matrix 
is  preceded  by  a  single  control  card  (Format  15)  which  indicates  the  number 
of  equipment  parameters  vectors  to  be  read  in  the  EN  matrix.    The  first  item 
read  in  an  equipment  parameters  vector  is  the  equipment  number.    The 
second  item  on  the  input  record  is  not  used  and  is  left  blank.    The  remaining 
thirteen  items  are  the  equipment  parameters  for  the  equipment  in  question. 
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The  equipment  parameters  are  unique  for  each  equipment  and  normally 
represent  such  parameters  as  power,  heat  transfer  coefficient,  etc.    The 
equipment  parameters  for  each  equipment  employed  in  the  steam  power 
system  are  listed  in  Figure  A9 . 

Equipment  Control  Parameters  Matrix  (ENC) .    Format:    5F15.5.    The 
ENC  matrix  is  preceded  by  a  single  card  (Format  15)  which  indicates  the 
number  of  rows  in  the  ENC  matrix.     Representative  data  to  be  used  in  the 
ENC  matrix  would  normally  include  allowable  differentials  across  equip- 
ments, costs,   etc.     Control  vectors  are  only  included  for  those  equipments 
which  employ  them.     If  no  control  vectors  are  employed,  the  control  card 
should  be  punched  zero  or  left  blank  and  input/output  should  proceed  to 
next  data  segment. 

Stream  Variables  Matrix  (SN) .    Format:    5FI5.5.    The  SN  matrix  is 
preceded  by  a  single  card  which  indicates  number  of  stream  vectors  to  be 
read  into  program.     If  control  variable  is  zero  or  blank,  the  input/output 
will  proceed  to  next  data  group.    The  stream  variable  matrix  contains  such 
variables  as  pressure,  temperature,  enthalpy,  total  steam  flow,  etc. 
This  location  labeling  is  consistent  for  each  subroutine.    The  SN  vector 
for  the  steam  power  system  is  structured  as  follows. 


STRM. 

FLAG 

PRESSURE 
(PSIA) 

TEMPERATURE 
(°F) 

ENTHALPY 
(BTU/LB) 

SPEC.  WT. 
(LB/FT3) 

VELOC. 
(FPS) 

TOTAL  FLOW 
(LB/HR) 

H20  FLOW 
(LB/HR) 

STEAM  FLOW 
(LB/HR) 

IN  PURITIES  FLOW 
(LB/HR) 

AIR  FLOW 
(LB/HR) 

—     ...    ....  _   . 

FUEL  FLOW 
(LB/HR) 

FLUE  GAS  FLOW 
(LB/HR) 

COOLANT  FLOW 
(LB/HR) 
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Stream  Control  Matrix  (SNC) .    Format:    5FI5.5.    A  single  card 
precedes  SNC  and  indicates  number  of  stream  control  variable  vectors  to 
be  read.     If  this  value  is  zero,   PACER  proceeds  to  next  data  group  in  data 
deck.    The  SNC  matrix  contains  control  parameters  such  as  viscosity, 
volume,  etc.    The  distinction  between  the  variables  specified  for  this 
matrix  and  that  provided  for  SN  are  arbitrary  and  lack  discrimination. 
Accordingly,  SNC  need  not  be  used  unless  required.    It  is  not  employed 
in  the  steam  power  cycle. 

Stream  Variables  Test  Vector  (DELS).    Format:    5FI5.5.    This  group 
is  preceded  by  a  single  card  (Format  15)  which  indicates  test  vectors  is 
to  be  read  from  DELS  matrix.    If  control  card  =  0,  no  test  vector  is  to  be 
read  and  PACER  skips  to  next  data  group.     If  control  card  value  ^  0  then 
PACER  will  read  one  test  vector  NSLMAX  long.    This  same  test  vector 
will  be  used  as  test  vector  for  determining  convergence  of  each  stream 
variable  for  all  streams. 

Stream  Control  Variables  Test  Vector  (DELSC) .    Format:    5FI5.5. 
Preceded  by  single  control  card  which  indicates  if  test  vector  is  employed. 
If  control  card  value  =  0  (Format  15)  then  no  control  test  vector  is  employed 
and  PACER  will  proceed  to  next  data  group.     If  control  card  value  7^  0, 
then  test  vector  is  employed  and  is  NSCLMX  long.     Only  one  test  vector 
is  employed  for  stream  control  variables  for  all  streams  employing  control 
variables.    This  test  vector  will  check  convergence  of  control  variables 
for  all  streams  employing  control  variables. 

Preferred  Stream  Numbers  (KPS) .    Format:    1015.     KPS  is  preceded  by 
a  control  card  (Format  15)  which  indicates  number  of  preferred  streams. 
If  control  card  value  equals  zero,  then  PACER  proceeds  to  next  data  group. 
If  control  card  7^  0,   PACER  reads  numbers  of  preferred  streams,  i.e.  those 
streams  which  will  be  assumed  known  first  in  a  trial  and  error  calculation 
to  determine  order  of  calculations.     Normally  these  should  be  recycle  or 
other  internal  streams  since  feed  and  product  streams  are  assumed  known 
in  trial  and  error  calculations. 

Extra  Equipment  Vectors  (AEN) .    Format:    5FI5.5.    AEN  is  preceded 
by  a  control  card.     If  control  card  value  (Format  15)  equals  zero,  PACER 
will  proceed  to  next  data  group.     If  control  value  ^  0,   PACER  reads  into 
AEN  extra  equipment  parameters  for  number  of  equipments  indicated  on 
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control  card.    This  matrix  is  only  used  if  insufficient  space  is  available 
in  EN  and  ENC.     Length  of  vector  is  arbitrarily  established  at  40. 

Order  of  Equipment  Calculation  (NELIST).    Format:    1015.     NELIST 
is  preceded  by  a  control  card  which  indicates  number  of  equipments  con- 
tained in  NELIST  matrix.     If  the  order  of  calculation  is  not  specified, 
the  control  card  value  (Format  1015)  equals  zero  and  PACER  proceeds  to 
next  data  group.     Each  card  consists  of  five  pairs  of  data.    The  first 
value  of  a  pair  represents  the  equipment  number.    The  second  value 
represents  the  calculational  mode  for  that  equipment.     If  mode  indicator 
equals  one,  the  equipment  will  be  calculated  directly.    If  the  mode 
indicator  equals  two,  the  equipment  will  be  calculated  as  an  equipment 
in  an  iterative  list  of  equipment.    Finally,  if  mode  indicator  equals  three, 
the  equipment  will  be  calculated  as  last  item  in  list  of  equipment  to  be 
calculated  in  an  iterative  mode. 

Selected  Streams  Vector  (KSPRNT) .    Format:    1615.     Contains  as  the 
first  item  the  number  of  stream  variable  vectors  to  be  printed  each  itera- 
tion.    If  it  is  not  desired  to  print  any  stream  variable  vectors  each  itera- 
tion the  card  should  be  left  blank.    The  stream  numbers  for  which  the 
stream  variables  are  to  be  printed  are  listed  after  the  first  item.    A  maxi- 
mum of  15  stream  variable  vectors  may  be  printed  each  iteration. 
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Appendix  B 

RESULTS  OF  TESTS 

Format  and  Interpretation  of  Results.    As  discussed  in  Chapter  V, 
six  simulation  tests  were  developed  to  test  the  validity  of  the  revisions 
made  to  the  original  PACER  program,  the  efficacy  of  PACER  as  a  simulation 
device  and  the  potential  usefulness  of  PACER  as  an  instructional  device 
in  a  chemical  engineering  curriculum.    The  results  of  those  simulation 
tests  are  presented  as  Figures  Bl.  -  B6. 

The  final  values  of  the  stream  variables  matrix  and  the  equipment 
parameters  matrix  are  presented  for  each  run.     In  addition  for  runs  01  and 
03  extensive  diagnostics  of  the  PACER  solution  process  are  presented 
for  analysis  by  the  user.    This  is  a  PACER  option.    The  stream  variables 
are  consistent  for  all  streams  in  the  stream  variables  matrix.    That  is, 
stream  variable  5  will  represent  enthalpy  for  all  streams.     However,  the 
enthalpy  for  all  streams  will  not  necessarily  be  the  same. 

For  the  steam  power  system  developed  for  this  project,  the  following 
stream  variable  arrangement  was  decided  on  and  should  be  used  as  a 
guide  in  interpreting  values  of  stream  variables  in  the  stream  variables 
matrix.     Note  that  there  are  five  columns  in  each  stream.    This  dimension 
may  be  parametically  varied  by  the  user  by  adjusting  dimensions  of  SN 
matrix. 
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0. 

0. 

16,0000 

0. 

5  o  o ,  n  o  n  n 

46?'.0fl00 
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STREAM  ViRU0l.ES  MATRIX  AT  RUN 

1    LOOP 

5 

1,0000 
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0. 

7,0000 

8. 

8, 

0.         169938,5684 
8,              8, 

0. 
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0.        640313?, OOOfl 
0,              8, 

0,            500,0000 
0.          2041,6575 

o.          a, 

70.0000 

8. 

0. 
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OUIPMENT  PARAMETERS  MATRIX 
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0. 

o. 

0. 
0, 
0. 

70,0000 

o, 

11 , 

93.0000 
0. 
0. 

99ft. 0000 
0. 
0. 

•i.0000 

0, 

0, 

0. 
0. 
0. 

1  ,0000 
0. 

o, 

0. 
0. 
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0. 

0. 

o, 
o. 
o, 

11 ,0000 

o. 

o , 

0, 
0, 

o. 

0,4800 

o. 
o. 

0, 
0. 
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STREAM  VARIABLES  MATRIX 
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1  •  FINAL 
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0, 

0, 

_  2,0000 
0, 

_#!___ 

5,0000 

6, 

8.4969 

4,0000 
0. 
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5,0000 
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8,4949 

6,0000 

o, 

n! 

■  7,0000 
0, 
0. 

8.000C 
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o. 

9,0000 
8, 
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8. 
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e, 
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8, 
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Of 
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v 
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o; 

or 
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0, 
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o, 
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0. 
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o, 
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0. 

0. 
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li 
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8, 
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0. 
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n, 
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0, 
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0. 
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0, 
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8 
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0. 

8, 
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3n 
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•  13  EOUtPMgHT  INTERCONNECTED  RV      36  STREAMS. 
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STREAM  VARIABLES  MATRIX 
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